An Angle Dependent SSC Model for Relativistic Jet Sources 
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ABSTRACT 

We report on the development of a numerical code to calculate the angle-dependent syn- 
chrotron + synchrotron self-Compton radiation from relativistic jet sources with partially ordered 
magnetic fields and anisotropic particle distributions. Using a multi-zone radiation transfer ap- 
proach, we can simulate magnetic-field configurations ranging from perfectly ordered (unidirec- 
tional) to randomly oriented (tangled). We demonstrate that synchrotron self-Compton model 
fits to the spectral energy distributions (SEDs) of extragalactic jet sources may be possible with 
a wide range of magnetic-field values, depending on their orientation with respect to the jet axis 
and the observer. This is illustrated with the example of a spectral fit to the SED of Mrk 421 
from multiwavelength observations in 2006, where acceptable fits arc possible with magnetic-field 
values varying within a range of an order of magnitude for different degrees of B-field alignment 
and orientation. 



Subject headings: radiation mechanisms: non-thermal 
galaxies — relativistic processes 

1. Introduction 

Blazars form one of the most energetically ex- 
treme classes of Active galactic nuclei (AGN). 
Blazars can be observed in all wavelengths, rang- 
ing from radio all the way up to 7-rays. Their 
spectral energy distribution (SED) is character- 
ized by two broad non-thermal components, one 
from radio through optical, UV, or even X-rays, 
and a high-energy component from X-rays to 7- 
rays. In addition to spanning across all observable 
frequencies, blazars are also highly variable across 
the electromagnetic spectrum, with timescales 
ranging down to just a few minutes at the highest 
energies. 

There are two fundamentally different ap- 
proaches to model the SEDs and variability of 
blazars, generally referred to as leptoni c and 
hadronic models (see, e.g., iBottcher 120071 for a 
review of blazar models). In the case of lep- 
tonic models, where leptons are the primary 
source of radiation, synchrotron, synchrotron self- 
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Compton (SSC), and external-Compton (EC) 
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present study is also on a leptonic model. In 
hadronic models, the low-energy SED component 
is still produced by synchrotron emission from 
relativistic electrons, while the high-energy com- 
ponent is dominated by the radiative output from 
ultrarelativistic protons, through photo-pion in- 
duced cascades and proton synchrotron emission 



(e.g.. Mannheim & Bicrmann 1992; Mucke fc Protheroe 
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200l[ iMiicke et all 120031 )7 One aspect common to 
all blazar models is a relativistic jet oriented at 
a small angle with respect to our line of sight, 
resulting in relativistic Doppler boosting and the 
shortening of observed variability time scales. 

Given computational limitations, the complex 
physical processes in relativistic jets can, realis- 
tically, only be evaluated with certain simplifying 
approximations. In order to facilitate analytical as 
well as numerical calculations, the two most com- 
mon approximations employed in blazar jet mod- 
els are to assume that the magnetic (B) field is 
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randomly oriented and tangled, and that the lep- 
ton momentum distribution is isotropic in the co- 
moving frame of the high-energy emission region. 
These two assumptions greatly simplify the evalu- 
ation of the synchrotron and Compton emission by 
eliminating various integrals over the interaction 
and scatterin g angles. However, there is increas- 
ing evidence (Attridge et al. 1999L Lvutikov et al. 
l2005l : lMarscher et al.ll2008l ) for a fairly well defined 
helical B-field structure within AGN jets. These 
observations also suggest a spine-sheath geometry 
for AGN jets. The differential velocity profiles 
within the jet is expected to create anisotropies 
in the particle distributions. It is therefore im- 
portant to explore jet models where we can not 
only simulate an ordered B-field, but also study 
the resulting radiation behaviour with anisotropic 
lepton distributions. 
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Fig. 1. — The Bi SO /B a ratio shows how the B- 
field estimate can vary depending on the elec- 
tron power-law index and whether one assumes 
an isotropic B-field or a specific pitch angle. 



1.1. B-Field estimates 

The standard approach to diagnosing the mag- 
netic field properties is via synchrotron polariza- 
tion. If the underlying distribution of emitting 
electrons is a power-law with power-law index p, 
the maximum degree of synchrotron polarization 
is given by: 
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where P± (y) and P» (v) are the synchrotron power 
per unit frequency in directions perpendicular and 
parallel to the projection of the magnetic field on 
the plane of the sky. Using equation [TJ we can see 
that for a power-law index of p = 3, the degree of 
polarization can be as high as 75%. 

It is therefore possible to estimate the magnetic 
field orientation based on polarization measure- 
ments, but an estimate of the field strength usu- 
ally requires the consideration of flux and spectral 
properties of the synchrotron emission. Further- 
more, polarization measurements are notoriously 
difficult (and even barely feasible at frequencies 
higher than optical), and may often not give re- 
alistic results due to Faraday rotation and depo- 
larization along the line of sight. In this work, we 
are interested in taking a complementary approach 
to estimating the magnetic field orientation where 
the difference in observed flux levels of the spec- 
trum can give an estimate of how the magnetic 
orientation may be changing. 



The principle behind this approach can be 
demonstrated when one compares the B-field es- 
timates based on a power-law distribution of elec- 
trons with an arbitrary power-law index, p , and 
pitch angle, a. The comparison of synchrotron 
emission coefficients for a power-law distribution 
of electrons with and without pitch-angle (a) de- 
pendence gives us a measure of how the estimated 
magnetic field strength can differ. The em ission 
coefficients can be found in lLongairl (|1994) and 
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Where k is the electron distribution power-law 
normalization. The above two expressions can be 
solved for the magnetic field to obtain: 
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This gives an estimate of how, for a given lumi- 
nosity, the B-field estimates can differ depending 
upon whether we assume an isotropic pitch angle 
approximation or a given pitch angle (which, in 
the case of relativistic electrons, is equal to the 
angle between the magnetic field and the line of 
sight). The above relation is only applicable in 
the optically thin regime. We can see in figure 
[T] that depending on the pitch angle assumption, 
and the electron distribution power-law index, the 
Biso/ B a fraction can range from (there is negli- 
gible synchrotron emission along an ordered mag- 
netic field) to ~ 1.25. Because the Compton emis- 
sivity is approximately isotropic for an isotropic 
distribution of electrons, the ratio Fssc/Fsy will 
change with the pitch angle. It is therefore im- 
portant to see how the overall synchrotron and 
synchrotron-self Compton spectra differ with well 
ordered magnetic fields. One point worth noting is 
that in our set-up the lower limit on B-field orien- 
tation is limited by the 5(fl sy — f2 e ) approximation 
(see section [2~2|) ; a magnetic field perfectly aligned 
with the observing direction will give zero output. 
However, in a more rigorous treatment, the lower 
limit on the minimum angle, a, for the B-field 
orientation will be determined by the relativis- 
tic beaming characteristic of synchrotron emission 
along an electron's direction of motion into a cone 
of opening angle a sy ~ 1/7, where 7 is the elec- 
tron Lorentz factor. In the case of optical frequen- 
cies and magnetic fields of B ~ 1 G, a m i n will 
be the order of 10~ 4 , while in the X-ray regime 
a-min ~ 10~ 5 . Using the relation in equation|4]it is 
possible to estimate the effects on the B-ficld esti- 
mates between this value of a m i„ and an isotropic 
magnetic field. For a power-law index of 3, a m j„ 
gives B iso /B a value of 0.0039. 

2. The Model 

The following section briefly outlines our model, 
including the synchrotron radiation and Compton 
scattering treatments followed as well as the nu- 
merical techniques used to implement them. 

2.1. Volume, B-field, and Distributions 

In our model the basic volume structure is a 
cubic cell. This allows the model to be modular 
and build an arbitrarily large volume with any de- 
sired anisotropies. Each cell contains a magnetic 
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Fig. 2. — Top: The electron and B-field interac- 
tion angles. Bottom: Cell/Jet geometry illustrat- 
ing how the photon spectrum can be observed as 
a function of various angles, 9 ce u and ce /;, with 
respect to the cell. The electron and photon dis- 
tributions are defined with respect to the cell. 

field plus electron and photon distributions. The 
magnetic field can have an arbitrary orientation 
and strength in each cell. This means that the 
overall volume can be modelled to contain a com- 
pletely uniform, partially anisotropic, or pseudo- 
random B-ficld. The purpose of the present work 
is to isolate the effects of the degree of order and 
orientation of the magnetic field on the emerging 
synchrotron emission. Therefore, we choose the 
simplest conceivable approach concerning the elec- 
tron distribution, and do not take electron cooling 
into account. This means, we only focus on static 
electron distributions which do not evolve due to 
energy losses. In future work we aim to include 
self-consistent cooling effects which would also al- 
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Fig. 3. — Volume/Jet geometry illustrating how a 
block of cells can be used to construct a larger vol- 
ume. Oview, the viewing angle with respect to the 
volume/jet is incorporated into Doppler boosting 
calculations. 



low us to probe how the pitch-angle dependence of 
the synchrotron cooling would give rise to different 
electron distributions in different cells, depending 
on the magnetic field set-up. 

The directional information in the electron and 
photon distributions and the B-field is with re- 
spect to the cell. In the case of electrons the dis- 
tribution is a function of energy and two angles 
with respect to the cell (see figured]). This gives 
us the ability to create anisotropies in the electron 
distributions as well by either having preferential 
direction for the electrons or by setting up the 
electron distributions differently in various cells. 
For both electrons and photons, the distributions' 
energy grids (Lorentz factor 7 for electrons and 
frequency v for the photons) are calculated using 
logarithmic binning. Therefore each distribution 
is modelled using a 3 dimensional array with the 
dimensions of [bins 7 /„ x binsg x bins^]. Figures [5] 
(right) illustrates how various angles with respect 
to the cell are defined. Angles 9 ce u and 4> ce u run 
from to 7r and to 2ir respectively. In figure[3]we 
can see how the overall volume can be constructed 
from individual cells. For a given viewing angle, 
the emission from the visible outer layer of cells is 
combined to produce an overall spectrum from an 
effectively larger volume. 



The simulation currently transfers, from one 
cell to another, only the photons. In order to 
achieve this, we need to calculate which of the six 
cubic faces a given photon direction will intersect. 
To calculate this, we assume that all the photons 
are produced in the center of the cell, and then 
trace photon paths in any given direction towards 
the nearest boundary. Even though our simula- 
tion considers a static situation, the transfer and 
radiative feedback between different cells requires 
an inherent time-dependence in the code. The 
time step for our radiation transfer approach is 
the light crossing time across a single cell, which 
is equal to the time it takes for the photons to 
travel from one cell to another. At the end of 
each time step and depending on the physical pro- 
cesses being modelled, the photon distribution is 
modified and passed to the appropriate neighbour. 
When being passed to a neighbour the entire pho- 
ton distribution is passed. Therefore at the end of 
a time-step each cell's (intrinsic) photon distribu- 
tion is emptied into six neighbouring cells, unless 
it is a boundary cell. The six incoming (transiting) 
photon distributions are stored until the start of 
the following time-step when they are combined to 
form a single intrinsic photon distribution again. 
The physical processes are then carried out on this 
single photon distribution. Synchrotron radiation 
is calculated first and the photons added to the in- 
trinsic distribution. Compton scattering is carried 
out after the synchrotron radiation. At this point 
we reach the end of a time-step and the process of 
transferring photon distributions to neighbouring 
cells begins again. The observed photon distribu- 
tion originates from the boundary cells. The pho- 
ton distributions emerging from visible faces of the 
boundary cells are combined to create a single ob- 
served photon distribution. This process of com- 
bining the photon distributions from the boundary 
cells in effect treats the whole multi-cell structure 
like a single cubic/cuboid structure. 

2.2. Synchrotron radiation 

Here we highlight the key points of the syn- 
chrotron radiation treatment that we follow. A 
m ore in-depth an alysis and details can be found 
in lLongairl (|l994l) . 



The synchrotron emissivity per electron, P v , is 
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given by: 
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where F(x) is given by: 
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x = v/v c , where v c is the critical frequency given 
by 37 2 eBsina/47rm e c. a, the pitch angle, is cal- 
culated using spherical trigonometry: 



cos(a) = cos 
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The synchrotron emission coefficient is given by: 

J Air Jl 
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The numerical Bessel function integration in equa- 
tion[5]can be time consuming. However, some fast 
routines to perfo rm this integration are given by 
Umstatterl (|l98lh which we modified for our pre- 
cision and computer language. 

In a full treatment of the synchrotron radiation 
the emitted photons are distributed within a solid 
angle (fi ~ I/7 2 ) about the pitch angle a. How- 
ever, for our purposes we assume the emitted pho- 
tons travel in the same direction as the emitting 
electrons. 

A detailed calculation of syn c hrotro n self- 
absorption can be found lLongairl (|1994l ). The 
absorption coefficient when recast in terms of elec- 
tron Lorentz factors, 7, instead of E, can be writ- 
ten as: 
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The photons produced via synchrotron radiation 
are added to the intrinsic photon distribution of 
the cell. The photons received from neighbouring 
cells are added to the intrinsic photon distribu- 
tion prior to calculating the synchrotron spectrum. 
Therefore the photons passing through any cell arc 
also synchrotron self-absorbed. The emission and 
absorption coefficients are used to calculate the 
total spectrum, 
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where r„ = \l is the optical depth and / is the size 
of the emission zone/cell. 

2.3. Compton scattering 
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Fig. 4. — Various angles involved in the Comp- 
ton scattering of photons with energy e to e s off 
electrons e 

In the limit 7 >> 1, and in the electron rest 
frame, the incident photon travels in nearly the 
opposite direction to the electron. This is due to 
photon aberration: 
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f3 e COS tp 

When /3 e — > 1 , cos 6' —> — 1 , we are in the head-on 
approximation regime, which we employ to greatly 
simplify Compton cross section calculations. That 
is, we can assume that the scattered photon solid 
angle, Q s , is well approximated by the electron 
solid angle f2 e . When the differential Compton 
(Klein- Nishina) cross section is i ntegrated over fi s , 
we get (jDermer fc Menonl [20091) : 
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where H is a Heaviside function, e' = 7e(l — 
flecosip), and the Compton kernel is given by: 

_i 2e, 



y + y 



(13) 
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and y = 1 — (e s /j). The Compton cross section 
can then be used in the emission coefficient for- 
mula to obtain the Comptonized spectrum. The 
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head-on approximation simplifies the emission co- 
efficient calculation by eliminating two integrals 
from the Compton emissivity treatment without 
the approximation. The following relation can be 
used to obtain the number of interacting photons: 
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The interacting photons are a combination of pho- 
tons originating from synchrotron radiation and 
the photons received from neighbouring cells. At 
the start of a time-step, the photon distributions 
received from the neighbouring cells are combined, 
while preserving the direction information, to form 
the intrinsic photon distribution. Synchrotron 
photons are also added to the intrinsic photon dis- 
tribution. The total photon distribution is then 
used in the Compton emissivity relation to obtain 
the Compton spectrum in the head-on approxima- 
tion, given by: 
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After the intrinsic photon distribution has been 
Compton scattered it is redistributed based on 
change in energy and direction. The redistributed 
photon distribution is then used to work out which 
neighbouring cells receive which proportion of the 
distribution. 

2.4. Overall spectrum 

The overall spectrum is obtained by combin- 
ing the synchrotron and Compton spectra. As 
it stands in our model, synchrotron emission is 
the only source of photons which are then Comp- 
ton scattered by the same population of electrons. 
The resulting photon spectrum is given as a func- 
tion of two cell angles ce u and 4> ce a (see figure 
Although not included in the present study, 
it is straight forward to include external Compton 
effects by adding the external photon field to the 
photon distributions. 



Once we have a spectrum, /„, the flux can be 
calculated using: 
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for an emission zone with an area A and luminosity 
distance dj,. For a viewing angle 0j the Doppler 
factor 8 if: is given by 



[T(1t/3(cos0,-)] 



(17) 



The 1= f' corresponds to either an approaching or 
a receding component of the jet. In the case of 
blazars the observed emission is strongly domi- 
nated by the approaching jet, boosted with the 
Doppler factor <5_. Also, any given frequency, is, 
in the emission region rest frame will be shifted by 
a factor of: 
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(18) 



We follow the I Wright) (120061) formulation to calcu- 
late the luminosity distance based on the redshift 
of an object. Photons emitted at an angle 0' h 
in the cell rest frame will appear at an angle p h 
due to angle aberration, which can be expressed 
as follows: 



cos p h 



cos0' ph + (3 



1 + P cos 
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(19) 



Results 



In order to study the effects of the B-field ori- 
entation on the synchrotron and synchrotron self- 
Compton spectra, we set up two scenarios. In the 
first set up the magnetic field is uni-directional in 
all the cells (27 in total) and in the second scenario 
the magnetic field is randomly oriented in each of 
the cells. These set-ups are likely to be the two ex- 
treme scenarios for a jet. Evidence points to the 
B-field being semi-ordered in AGN jets; for exam- 
ple, helical (Attridge et al. 19991 : Lvutikov et al. 
20051 : iMarscher et al.ll2008l ) . In all the presented 



cases, the electron distribution is a power-law and 
distributed uniformly over the angles and <f>. 

Figure [5] shows synchrotron self-Compton spec- 
tra with different B-field configurations. Various 
simulation parameters can be found in tabic [1] 
The spectra are for a fixed viewing angle. Chang- 
ing the B-field orientation has a significant impact 
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Fig. 5. — Synchrotron self-Compton spectra at 
a viewing angle of 0.014° (0.57° in jet frame). 
The spectra shown results of different B-field con- 
figurations, including when it is randomly ori- 
ented in different cells. The plot shows Ob val- 
ues (tpB — 90°, unless randomly oriented). The 
red dot-dot-dashed curve shows the correspond- 
ing calculation based on the angle averaged emis- 
sivity and spherica l geom etry, using the code of 
Bottcher fc Chiang! (|2002h . See figure [2] for details 
on various angles. 

on the observed spectrum. There is a large differ- 
ence in flux values depending on the B-field orien- 
tation with respect to the observing direction. The 
minimum flux levels are observed at an orientation 
along the viewing angle while the maximum flux 
levels are observed when the B-field is perpendic- 
ular to the line of sight. The figure also demon- 
strates the fact that only the synchrotron spec- 
trum component is heavily affected by the B-field 
orientation. The Compton scattered component 
of the spectrum is almost independent of the B- 
field orientation. The main reason for this is that 
the photon distribution anisotropies introduced by 
the B-ficld orientation are lost when scattering off 
an isotropic distribution of electrons. The small 
variations that remain are due to the 1 — /3 cos ip 
factor present in Compton emissivity calculations. 
The line of sight and photon anisotropies therefore 
affect the extent to which the Compton spectrum 
is boosted. Additional anisotropies are introduced 
by the discretization of the photon and the elec- 
tron angular distributions. 

Figure [5] also s hows a comparison with a nother 
SSC calculation ([Bottcher fc Chiang! liooH ) which 



Fig. 6. — Synchrotron self-Compton spectra at a 
various viewing angles (jet frame values). The B- 
field orientation: 9b — 45°, 4>b — 90°. See figure 
[2] for details on various angles. 



assumes a randomly oriented magnetic field. This 
simulation is set up with identical parameters to 
the ones outlined in table [TJ except it uses a spher- 
ical volume instead of a cube; the total volumes 
are identical, therefore the sphere has a radius of 
1.86 x 10 14 cm. We can see that the synchrotron 
components are in good agree ment. However, the 
Comp ton component in the Bottcher fc Chiand 
(|2002l ) calculation is much higher. The inverse 
Compton to synchrotron ratio Fjc / F syn differs by 
a factor ^0.68 between the two calculations. This 
is most likely due to differing geometries. For a 
sphere, the average photon escape time is 3i?/4c, 
whereas in our cubic set up, due to the way radia- 
tion transport between cells is treated (see Section 
12. ip . a photon takes, on average, 2d ce a/c to escape 
the region (d C M is the width of an individual cell) . 



Since the flux ratio Fgsc I 'Fsyn — u syn /uB and the 
volume-averaged radiation energy density u syn is 
proportional to the photon escape time scale, the 
longer photon escape time scale in the spherical 
geometry results in a larger SSC flux. 

Figures [6] shows the effects of the viewing angle 
on the spectrum. It shows SEDs for a fixed az- 
imuthal angle, but different viewing angles with a 
single B-field configuration. As before, there are 
normalization differences between the SEDs when 
comparing uniform and randomly oriented B-field, 
but the main factor in this case is the variation in 
the Dopplcr boosting due different viewing angles. 
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Parameter 




Cells 


27 


Cell size 


1 x 10 14 cm 


z 


0.031 


r 


20.0 




9 ergs cm -3 


p 


3.5 


imin 


5 x 10 3 


j max 


5 x 10 5 


B 


0.9 G 



Table 1: The simulation parameters for the results 
show in figures [5] and [6j Key to symbols: z = 
redshift, T = Jet Lorentz factor, u e = electrons 
energy density, p electrons distribution power law 
index, r y m in,max = electrons distribution minimum 
and maximum energy, B = B-field strength 



3.1. Markarian 421 
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Fig. 7. — Synchrotron self-Compton fit for 
Markarian 421. The spectral fits are at various 
B-field orientations. The data are from XMM- 
Newton OM, XMM-Newton EPIC, and VERI- 
TAS. See table [2] for the fit parameters. 

Markarian 421 was the first extragalactic source 
to be detected in TeV energies, hence making it 
an extensively studied source. We present spec- 
tral fits t o XMM-Newton and VERITAS data pre- 
sented in lAcciari et al. ( 20091 ). We refer the reader 
to the above paper for the details on data reduc- 
tion. 

Figures [71 [SI and H show fits to Mrk 421 data 
using our code. Our aim in the present paper 
is to explore the impact of the B-field orienta- 
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Fig. 8. — Synchrotron self-Compton fit for 
Markarian 421. The spectral fits are at various 
B-field orientations. The data are from XMM- 
Newton OM, XMM-Newton EPIC, and VERI- 
TAS. See table [2] for the fit parameters. 



tion on the fit parameters, especially the B-field 
estimates. In figure [7] we can see that the syn- 
chrotron component is best fit for a particular B- 
field orientation (see table [2] for fit parameters). 
The gamma-ray data, however, is fit well with 
all the orientations. This is due the fact that 
the self-Compton spectrum is not affected much 
by the B-field orientation (see the discussion in 
the previous section). There is a significant dif- 
ference in the synchrotron peak when comparing 
various B-field orientations. The spectra in fig- 
ure [7] show that at a good fit is achieved when 
the B-field is oriented at 85° with a strength of 
0.18 G. However, the fit shown is not unique. We 
can see in figures |5J and [5] that for identical pa- 
rameters, except the B-field strength, the best-fit 
B-field orientation is very different. In one case a 
pseudo-random B-field provides the best fit with 
B = 0.22 G, while for B = 0.25 G, a B-field orien- 
tation of 45° provides the best fit. Therefore good 
fits can be achieved with different B-field strengths 
at different orientations. The main point here is 
the fact that it is possible to over- or underesti- 
mate the B-field strength when assuming it to be 
randomly oriented. In the cases presented here, 
the B-field strength ranges from 0.18 G to 0.25 G 
for very similar fits to the data, but with different 
magnetic-field orientations. Therefore it is pos- 
sible to overestimate the magnetic field strength 
by at least 30% if a particular B-field orientation, 
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Fig. 9. — Synchrotron self-Compton fit for 
Markarian 421. The spectral fits are at various 
B-field orientations. The data are from XMM- 
Newton OM, XMM-Newton EPIC, and VERI- 
TAS. See table [2] for the fit parameters. 

whether uniform or tangled, is assumed. If the 
B-hcld were pointed closely aligned with the line 
of sight, much higher B-field values will be neces- 
sary to obtain similar fits (see discussion in section 
II. ip . We also note that the bulk Lorentz factor 
used in the fits are lower than the values obtained 
by some authors fo r fitting Mrk 421 data (e.g. see 
Aleksic et al.ll2011 ) . The value of the bulk Lorentz 



factor values used in our fits is likely to be on 
the lower end of the limits impose d by pair opac- 
ity arguments (jCelotti et al.lll998l ). However, the 
main point of this paper is not the determination 
of actual best-fit values for Mrk 421 (which would 
not be realistic due to our neglect of cooling ef- 
fects anyway) , but to demonstrate the orientation- 
dependent magnetic-field degeneracy in the course 
of blazar SED fitting. 

4. Conclusions 

In this paper we have presented first results 
from a new relativistic jet radiation transfer code 
that we are currently developing. Here we take the 
full angular dependence into account when mod- 
elling synchrotron and synchrotron self-Compton 
processes. We are able to model the B-field at ar- 
bitrary orientations and study its impact on the 
resulting spectra. 

We have seen that the B-field orientation plays 
an important role on the normalization of the 
synchrotron spectrum. Using fits to Markarian 
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Table 2: Fit parameters for the Markarian 421 data 
shown in figures [7] and HI Key to symbols: z = 
redshift, T = Jet Lorentz factor, u e — electrons 
energy density, p electrons distribution power law 
index, ^min.max = electrons distribution minimum 
and maximum energy, B — B-field strength 



421 data, we have shown how the B-field orien- 
tation can mislead into over/under-estimating its 
strength. Any future work should therefore be 
mindful of the fact that the underlying assumption 
about the B-field orientation will play a consider- 
able role in the errors associated with the magnetic 
field strength estimates. 
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